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Abstract. 

In this paper, we review the current theory of very low mass stars 
model atmospheres including the coolest known M dwarfs, M subdwarfs, 
and brown dwarfs, i.e. T^s < 5,000K and -2.0 < [M/H] < +0.0. We 
discuss ongoing efforts to incorporate molecular and grain opacities in 
cool stellar spectra, as well as the latest progress in deriving the effective 
temperature scale of M dwarfs. We also present the latest results of the 
models related to the search for brown dwarfs. 



1. Very Low Mass Star models and the Tofj Scale 

Very Low Mass stars (VLMs) with masses from about 0.3 Mq to the hydrogen 
burning minimum mass (0.075 Mq, Baraffe et al. 1995) and young substellar 
brown dwarfs share similar atmospheric properties. Most of their photospheric 
hydrogen is locked in H2 and most of the carbon in CO, with the excess oxygen 
forming important molecular absorbers such as TiO, VO, FeH and H2O. They 
are subject to an efficient convective mixing often reaching the uppermost layers 
of their photosphere. Their energy distribution is governed by the millions of ab- 
sorption lines of TiO and VO in the optical, and H2O in the infrared, which leave 
no window of true continuum. But as brown dwarfs cool with age, they begin 
to differentiate themselves with the formation of methane (CH4) in the infrared 
(Tsuji et al. 1995) at the expense of CO which bands begin to weaken in their 
spectra (Allard et al. 1996). Across the stellar-to-substellar boundary, clouds 
of e.g. corundum (AI2O3), perovskite (CaTiOa), iron, enstatite (MgSiOs), and 
forsterite (Mg2Si04) may form, depleting the oxygen compounds and heavy el- 
ements and profoundly modifying the thermal structure and opacity of their 
photosphere (Sharp & Huebner 1990, Burrows et al. 1993, Fegley & Loggers 
1996, Tsuji et al. 1996ab). 

Because these processes also occur in the stellar regime where a greater 
census of cool dwarfs is currently available for study, a proper quantitative un- 
derstanding of VLM stars near the hydrogen burning limit is a prerequisite to an 
understanding of the spectroscopic properties and parameters of brown dwarfs 
and jovian-type planets. Model atmospheres have been constructed by several 
investigators over recent years with the primary goals of: 



1. Determining the effective temperature scale of M dwarf stars down to the 
substellar regime. 

2. Identifying spectroscopic signatures of substellarity i.e. gravity indicators 
for young brown dwarfs, and spectral features distinctive of cooler evolved 
brown dwarfs. 

3. Providing non-grey surface boundary to evolution calculations of VLMs 
and brown dwarfs leading to more consistent stellar models, accurate mass- 
luminosity relations and cooling tracks for these objects. 

The computation of VLMs and brown dwarf model atmospheres requires 
a careful treatment of the convective mixing and the molecular opacities. The 
convection must currently be handled using the mixing length formalism while 
a variety of approximations have been used to handle the millions of molecular 
and atomic transitions that define the spectral distributions of VLMs and brown 
dwarfs. The most accurate of these methods is the so-called opacity sampling 
(OS) technique which consists in adding the contribution all transitions absorb- 
ing within a selected interval around each point of a pre-determined wavelength 
grid (typically ~ 22000 points from 0.001 to 100 /xm). When the detail of the list 
of transitions is lacking for a molecule as is the case for the important absorber 
VO, the Just Overlapping Line Approximation (JOLA) offers an alternative by 
approximating the band structure based on only a few molecular rotational con- 
stants. The straight-mean (SM) and K-coefficients techniques, which consist 
in averaging the opacities over fixed wavelength intervals chosen smaller than 
the resolution of typical observations, have also been used in modeling late-type 
dwarf atmospheres. Their main advantage is to save computing time during the 
calculation of the models, often at the expense of an accurate spectral resolution. 
The list of recent model atmospheres and the opacity technique they mostly rely 
upon is given in table |l|. 



Table 1. Relevant Model Atmospheres 



Authors 


Grid 


Teff range (K) 


Main Opacity Treatment 


Kurucz 1992 


Atlasl2 


3500 - 




OS 


AUard 1990 


Base 


2000 - 


3750 


SM+JOLA 


Saumon et al. 1994 


zero-metallicity 


1000 - 


5000 


OS 


Tsuji et al. 1995 


grainless 


1000 - 


2800 


JOLA 


Brett 1995 


MARCS 


2400 - 


4000 


OS 


AUard & Hauschildt 1995 


Extended Base 


1500 - 


4500 


SM 


Tsuji et al. 1996 


dusty 


1000 - 


2800 


JOLA+Grains 


AUard et al. 1996 


Next Gen 


900 - 


9000 


OS 


AUard et al. 1997b 


NextGen-dusty 


900 - 


3000 


OS+Grains 


Marley et al. 1996 






1000 


K-coefficients 



Because they mask emergent photospheric fluxes that would otherwise es- 
cape between absorption lines, the JOLA and SM approximations generally led 
to an excessive entrapment of heat in the atmosphere which yields systemati- 
cally hotter model structures, and higher effective temperature (Teg) estimates 
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Figure 1. Current model-dependent effective temperature scales for 
cool stars down to the hydrogen burning limit. Triangles feature results 
from spectral synthesis of selected stars from the works of Kirkpatrick 
et al. (1993) and Leggett et al. (1996) as indicated. The new generation 
of OS models by Brett (1995b) and Allard et al. (1996), as interpolated 
onto theoretical isochrones by Chabrier et al. (1996), reproduce closely 
the independently-determined positions of the eclipsing M dwarf binary 
system CM Dra and YY Gem, and the empirical Tcfr scale of Jones et 
al. (1994). 



for individual stars. Allard et al. (1997) have reviewed in detail the results of 
brown dwarfs and VLM model atmosphere calculations with respect to the effec- 
tive temperature scale of M dwarfs. We reproduce in Figure |l|the T^s — — I) 
relation of Allard et al. (1997) for the models listed in Table]^. 

Two double-line spectroscopic and eclipsing M dwarf binary systems, CM 
Draconis and YY Geminorum, offer some guidance in the sub-solar mass regime 
and are reported in Figure |^ according to Habets & Heintze (1981). The use of 
an OS treatment of the main molecular opacities, in particular for TiO, appears 
to yield a break-through in the agreement of Tofj scales with these two M dwarfs 
binary system. The NextGen and MARCS models yield effective temperatures 
that are coincidentally in good agreement with those derived empirically from 
the H2O opacity profile by Jones et al. (1994) Q Note, however, that the Atlasl2 
OS models suffer from an inaccurate TiO absorption profile and a complete lack 
of H2O opacities, and are therefore clearly inadequate in the regime of VLM 
stars (i.e. below T^s ~ 4500 K) where molecular opacities dominate the stellar 
spectra and atmospheric structures. 

Some uncertainties on the metallicity of the CM Draconis system may soon 
disqualify the latter as a member of the disk main sequence (Viti et al. 1997). 
This stresses the importance of finding other low-mass eclipsing binary systems 
in the disk. These are hopefully soon to be provided by the 2MASS and DENIS 
surveys (see D. Kirkpatrick and X. Delfosse elsewhere in this volume). Much 
uncertainty remains, therefore, at the lowermost portion of the main sequence. 
The inclusion of grain formation (as discussed below) and more complete opaci- 
ties of TiO promise a better understanding of the stars and brown dwarfs in the 
vicinity of the hydrogen burning limit (the location of which is roughly indicated 
in Figure |l| by the termination point of the Allard et al. 1996 model sequence), 
but still remain to be ascertained. 



2. The Infrared Colors of Brow^n Dwarfs 



The DENIS and 2MASS infrared sky surveys will soon deliver large data bases 
of red dwarfs, brown dwarfs and perhaps extrasolar planets, which will necessi- 



^Note that a comparison to observed spectra reveals uncertainties of the order of 0.2 to 0.5 
mag on the pubhshed / magnitudes of the latest-type M dwarfs G1406, VBIO and LHS2924 
analyzed by Jones et al. (1994) and reported on Figure |l 
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Figure 2. The observed infrared spectral distribution of the dM8e 
star VBIO as obtained at UKIRT by Jones et al. (1994) (bold full hue) 
is compared to model spectra obtained using (from bottom to top): (i) 
the SM laboratory opacity profile of Ludwig (1971), (ii) the 20 million 
line list by J0rgensen et al. (1994), (iii) the preliminary ab initio line 
list of 6.2 million transitions by Miller & Tennyson (1994), and (iv) the 
latest ab initio list of 300 million lines by Partridge & Schwenke (1997). 
The models (shown as dotted lines) are all fully converged and normal- 
ized to the observation at 1.2 /xm. Their parameters were determined 
from a fit to the optical stellar spectra (not shown) and are nearly the 
same in all four cases. Note that all 300 million lines of the Partridge 
&; Schwenke list have been included in the model construction! 

Figure 3. The most recent models of late type dwarfs are compared 
to the photometric observations of field stars and brown dwarfs, and to 
Pleiades objects including the brown dwarfs PP115, Teidel and Calar3. 
Unresolved binarity is reflected in this diagram by a red excess in J—K. 
The red dwarfs newly discovered by DENIS (see X. Delfosse elsewhere 
in this volume) are also shown, although their photometry is still very 
uncertain at this point. The field brown dwarf Gliese 229B is off the 
scale to the blue in J — K due to strong CH4 absorption in the K 
bandpass. This diagram offers excellent diagnostics to identify brown 
dwarf candidates of the field (very red in either J — K 01 I — J) or of 
the halo (very blue in both I — J and J — K). 

tate the best possible theoretical foundation. A proper understanding of their 
colors is essential in the search for brown dwarfs. Brown dwarfs and giant plan- 
ets emit over 65 radiation in the infrared (> 1.0/xm). Yet the main difficulties 
met by VLMs and brown dwarf modelers in recent years has been to reproduce 
adequately the infrared (1.4 to 2.5;um) spectral distribution of dwarfs with spec- 
tral types later than about M6. All models listed in the central part of table |l| 
underestimate the emergent flux, most as much as 0.5 mag at the K bandpass, 
despite the different opacity sources used by the authors. Allard et al. (1994, 
and subsequent publications) have explored water vapor opacity data from vari- 
ous sources. Figure Q summarizes these results. Clearly, the water vapor opacity 
profile is quite uncertain and has varied with the degree of completeness and the 
assumptions used in the construction of the molecular model and its potential 
surface. The most recent and complete line list of Partridge & Schwenke suc- 
ceeds for the first time in reproducing the 1.6 //m opacity minimum, in the H 
bandpass, well enough for the atomic Nal resonance line to finally emerge in 
the synthetic spectrum, matching the observed feature. However, it fails to pro- 
vide an improvement in the K bandpass where the less complete list of Miller 
&; Tennyson still yield the best match of the models to the observed spectra. 
The NextGen models of Allard et al. (1996) are computed using the Miller &: 
Tennyson line list and are the only models to provide a match to the infrared 
colors of VLMs. This is shown in Figure ^ where the complete series of NextGen 
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models — as interpolated on the Baraffe et al. (1997) isochrones for 10 Gyrs 
and 120 Myrs and ranging from metallicities of [M/H]= —2.0 to 0.0 — are com- 
pared to the photometric field dwarfs' samples of Leggett (1992), Tinney et al. 
(1993), and Kirkpatrick et al. (1995). Other models series including those of 
Brett (1995) and the Extended grid of Allard & Hauschildt (1995, not shown) 
are distinctively bluer than the observed sequence, while the 10 Gyrs NextGen 
models of solar metallicity follow closely the empirical sequence^ of Kirkpatrick 
& McCarthy (1994) until spectral types of M6 (i.e. J - K ^ 0.85). Beyond 
this point, all models fail to reproduce the bottom of the main sequence into 
the brown dwarf regime as defined by G1406, VBIO, BRI0021 and GD165B. The 
models catch up only at the much lower Tefj of the evolved brown dwarf Gliese 
229B, i.e. 900-1000 K (Allard et al. 1996, Marley et al. 1996). 

The cause of the model discrepancies at the stellar-to-brown dwarf bound- 
ary can only be one that affects the cooler models for Gliese 229B in a far lesser 
obvious extent. Since the infrared spectral distribution is sensitive to the mixing 
length, yet without allowing for an improved fit of VLMs spectra, Brett (1995a) 
suggested that the problem lie in the inadequacy of the mixing length formalism 
for treating the convective transport in an optically thin photospheric medium. 
These concerns may also be augmented by uncertainties about the extent of 
the overshooting phenomenon in VLMs (see F. D'Antona elsewhere in this vol- 
ume). The convection zone recedes gradually below the photosphere as the mass 
(and Tgfj) decreases along the isochrones. This implies that the lithium test of 
substellarity (Rebolo et al. 1992) — which relies on the assumption that the 
brown dwarf is still fully convective and mixing lithium from its core to its pho- 
tospheric layers after 10^ yrs of age — is inapplicable for objects cooler than 
Teff < 2200 K. The presence of lithium in the spectra of a late-type (>M10) 
field dwarfs, if detected, could only reflect their initial abundances and not their 
substellar nature. The shrinking of the convection zone also allows a very good 
agreement between the models of Marley et al. (which includes adiabatic con- 
vection only for the optically thick layers of the atmosphere) and the models of 
Allard et al. (1996) (based on a more careful treatment of convection with the 
mixing length formalism) for the brown dwarf Gliese 229B (see Figure 5 of Al- 
lard et al. 1997). Yet the maximum radial extent of the convection zone occurs 
at around Teg = 3000 K, while the discrepancy with the infrared observations 
increases steadily towards the bottom of the main sequence. 

A more promising answer to the so called "infrared problem" may rather 
be found in the formation of dust grains in the very cool (typically Tiayer ~ 
Tcfj — 1000 K) upper layers of red and brown dwarf's atmospheres. Tsuji et al. 
(1996a) proposed, based on their results of including the effects of the forma- 
tion and opacities of three grain species (AI2O3, Fe, and MgSiOs) in their new 
"dusty" models, that the greenhouse heating of grain opacities, the resulting 
enhanced H2O dissociation, and the infrared flux redistribution can explain the 
infrared spectra of cool M dwarfs. The formation of perovskite dust grains at the 
expense of TiO may also explain the observed saturation (and disappearance in 



^Note that this sequence was defined by stars selected from their optical spectroscopic properties. 
The somewhat erratic aspect of the sequence in this infrared diagram reflects uncertainties in 
the photometry and perhaps in the age of the selected stars. 
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GD165B and Gliese 229B) of the TiO bands in the optical spectra of late-type 
red dwarfs (see also Jones & Tsuji elsewhere in this volume). The implications 
of this result is far reaching. Field brown dwarf candidates such as BRI0021 
and GD165B can be far cooler and less massive than previously suspected (see 
e.g. the NextGen-dusty model predictions in Figure|l]). If grains also form in 
the young Pleiades brown dwarfs PP115, Teidel and CalarS (Tcfj ~ 3000, 2800, 
and 2700 K respectively), lithium abundances derived from grainless models and 
synthetic spectra such as those of Pavlenko et al. (1995, see also elsewhere in 
this volume) may be overestimated, and the masses attributed to these objects 
possibly underestimated. Evolution models of brown dwarfs, which are sensitive 
to the treatment of the atmospheres (Baraffe et al. 1995, Chabrier et al. 1996), 
and their predicted Mass-lithium abundance and Mass-Luminosity relations may 
also be affected. 

And indeed, the temperatures and pressure conditions of the outer layers of 
red dwarfs are propice to the formation of dust grains as demonstrated years ago 
by Sharp &: Huebner (1990) and Burrows et al. (1993). However it was not clear 
at the time if the inward radiation of an active chromosphere, or the efficient 
convective mixing from the interior, would heat up these upper photospheric 
layers and disable grain formation. Another concern is that, under the gravities 
prevailing in M dwarfs, gravitational settling may occur that would eliminate 
large grains and their opacities from the photospheres over relatively short time 
scales. These possibilities still need to be thoroughly investigated, but clearly, 
grain formation is a process that occurs in M dwarf and brown dwarf model 
atmosphere and it must included in such calculations. 

In order to investigate which grains may form in the upper layers of M 
dwarfs, Allard et al. (1997b, in preparation) have modified the equation of 
states used in the NextGen models to include the detailed calculation of some 
1000 liquids and crystals, using the free Gibbs energies compiled by Sharp & 
Huebner. Their results showed that, besides the three species considered by 
Tsuji et al., the M dwarfs atmosphere were rich in condensates with Zr02, 
CasAlaSiOr, CaaMgSiOT, MgAl204, TiaOs, Ti407, CaTiOg, and CaSiOs show- 
ing up in models as hot as T^s = 2700 — 3000 K (i.e dM8-dM6)! The preliminary 
NextGen-dusty models have been computed using a continuous distribution of 
ellipsoid shapes and interstellar grain sizes (between 0.025 and 0.25 /xm) for the 
treatment of the opacities of the AI2O3, Fe, MgSiOs, and Mg2Si04 dust grains 
(see Allard & Alexander elsewhere in this volume for computational details). 
This contrast with the assumption of spherical grains with 0.1 fim diameters in 
the dusty models Tsuji et al. Both model sets are shown in Figures |l] and ^ 
As can be seen, the dusty models of Tsuji et al. provide the correct tendency 
of the coolest models to get rapidly very red (as much as J — K = 1.65 for 
GD165B) with decreasing mass for a relatively fixed I — J color. Those models 
are however systematically too red in I — J by as much as 1 mag and do not 
reproduce even the most massive M dwarfs while over-predicting the effects of 
grains in Gliese 229B type brown dwarfs (Tsuji et al., 1996b), a problem which 
must be related to the use of the JOLA treatment of molecular opacities in these 
models (see section ^ above). The NextGen-dusty models, on the other hand, 
show the onset of grain formation effects by a progressive deviation from the 
grainless NextGen models for J — K > 0.85, bringing an improved agreement 
with the observed sequence in the region where the grainless NextGen models 
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deviate. Of course, much remains to be improved in the computation of models 
with dust grains. The size distribution of various grain species, in particular 
those of the perovskite CaTiOs which is responsible for the depletion of TiO 
from the optical spectra of late-type dwarfs (eg. GD165B, see D. Kirkpatrick 
elsewhere in this volume) and of corundum (AI2O3) which accounts for most of 
the grain opacities in current models, is unknown for the conditions prevailing 
in M dwarfs atmospheres. It is conceivable that grains form more efficiently 
in M dwarfs atmospheres than in the interstellar medium and therefore their 
opacities are larger than considered in the NextGen-dusty models. We may as 
well be missing a number of important contributors (e.g. Zr02) to the total 
grain opacities in the models. Further investigations including time dependent 
grain growth analysis will be required to determine the true contribution of dust 
grains to the infrared colors of red and brown dwarfs. 

In the meanwhile, diagrams like that of Figure |3| may help in distinguishing 
interesting brown dwarfs candidates from large data banks of detected objects, 
and in obtaining an appreciation of the spectral sensitivity needed to detect new 
brown dwarfs. Models (Tsuji et al. 1995, Allard et al. 1996, Marley et al. 1996) 
and observations of Gliese 229B (see B. Oppenheimer elsewhere in this volume) 
have shown that methane bands at 1.7, 2.4 and 3.3 iim appear in the spectra 
of cool evolved brown dwarfs, and cause their J — K colors to get progressively 
bluer with decreasing mass and as they cool over time. Yet their I — J colors re- 
main very red which allows to distinguish them from hotter low-mass stars, red 
shifted galaxies, red giant stars, and even from low metallicity brown dwarfs that 
are also blue due to pressure-induced II2 opacities in the H—to—K bandpasses. 
Fortunately, grain formation and uncertainties in molecular opacities are far 
reduced under low metallicity conditions ([M/H]< —0.5). Therefore, model at- 
mospheres of metal-poor subdwarf stars and halo brown dwarfs are more reliable 
than their metal-rich counterparts at this point. This has been nicely demon- 
strated by Baraffe et al. (1997) who reproduced closely the main sequences of 
globular clusters ranging in metallicities from [M/H]= —2.0 to —1.0, as well 
as the sequence of the Monet et al. (1992) halo subdwarfs in color-magnitude 
diagrams (see G. Chabrier elsewhere in this volume). The colors of halo brown 
dwarfs as predicted by the NextGen models are therefore of quantitative quality 
await confrontation with the infrared colors of metal-poor subdwarfs from e.g. 
the Luyten catalog and the US Naval Observatory surveys. The sensitivity of the 
I — J index to the chemical composition of the atmosphere (clearly illustrated 
by the NextGen model grid) allows to distinguish brown dwarf populations in- 
dependently of an accurate knowledge of the parallaxes or distances involved. 
Even young brown dwarfs of lower gravity appear to form a distinct sequence at 
bluer I — J (and redder J — K) values then that of their older field star counter- 
parts as also evident from a comparison of the 10 Gyrs and 120 Myrs NextGen 
models. This gravity effect, and perhaps enhanced grain formation, may explain 
the scatter of spectroscopic properties observed among field dwarfs at the bot- 
tom of the main sequence (Kirkpatrick, this volume), as well as the systematic 
differences between Pleiades brown dwarfs and older field stars of same spectral 
type (i.e. same VO band strengths) noted by Martin et al. (1996). 

Gravity effects have also been found to affect the infrared spectra of cool 
evolved brown dwarfs such as Gliese 229B: Allard et al. (1996) reported a strong 
response of the 2.2 |Um opacity minimum to gravity changes which allowed to 
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Figure 4. Predicted absolute fluxes of brown dwarfs at 50 pc as com- 
pared to the sensitivity of ground and space-based platforms which will 
be or are currently applied to the search for brown dwarfs and extra- 
solar planets. The latter are values reported for the 5 a detection of 
a point source in 1 hr of integration, except for the three NICMOS 
cameras where the integration is limited to 40 minutes (Saumon et 
al. 1996). Models of both Allard et al. (1996) (fuh) and Marley 
et al. (1996) (dotted) are shown which simulate (i) a brown dwarf 
near the hydrogen burning limit (topmost spectrum: T^fi = 2000K), 
(ii) an evolved brown dwarf similar to Gliese 229B (central spectra: 
Tgfr = 900K and 960K) , and (iii) a brown dwarf closer to the deuterium 
burning limit (lowermost spectrum: Tefr = 500K). The corresponding 
black-body (dashed) are also shown for comparison. 

restrain the mass of the brown dwarf ^. The general spectral distributions of 
cool evolved brown dwarfs are well reproduced by current models despite the 
difference in their respective modeling techniques, and despite the uncertain- 
ties tied to grain formation and incomplete opacity data base for methane and 
ammonia. The models of Allard et al. (1996) and Marley et al. (1996) are 
compared in Figure ^ which also summarizes the predicted absolute fluxes that 
free-floating brown dwarfs would have at a distance of 50 pc. As can be seen, 
there is no clear cut distinction between brown dwarfs and planets; molecular 
bands most gradually form (dust, H2O, CH4 and NH3) and recede (TiO, VO, 
FeH, and CO) from the stellar to the planetary regime as the atmospheres get 
cooler. They remain very bright in the UK region, and become gradually red- 
der in the near-infrared I to J bandpasses, which allows their detection from 
ground-based facilities. Layers of dust clouds in their upper atmospheres may 
increase the albedo of extrasolar planets and cool brown dwarfs sufficiently to 
reflect the light of a close-by parent star, becoming therefore resolvable in the 
optical where the clouds are densest but the parent star is however brightest. 
The peak of their intrinsic spectral energy distribution is located at 4.5 ^m. 
At 5 fim, the hotter (younger or more massive) brown dwarfs and stars show 
strong CO bands which cause their flux to drop by nearly 0.5 dex relative to 
that at 4.5 fim. And between 4.5 and 10 //m, opacities of CH4 (and H2O in 
the hotter brown dwarfs) cause the flux to drop by 0.5 to 1.0 dex. Searches 
in the 4.5-5 fiui region should therefore offer excellent possibilities of resolving 
brown dwarfs and EGPs in close binary systems, and to find free-floating brown 
dwarfs if space-telescope time allocations allow. The detection limits of current 
and planned ground-based and space-based telescopes Saumon et al. (1996) are 
also indicated in Figure ^ which show that brown dwarfs within 50 pc would be 
easily detected by SIRTF in the 4.5-5.0 fim region. The drop in sensitivity of 
the various instruments redwards of 10 nm implies, however, that brown dwarfs 



■^Only within the error on the flux cahbration of the observed spectra which are unfortunately 
large for this object. 
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and planets cooler than Gliese 299B have little chance to be detected in those 
redder bandpasses. 

3. Conclusions 

In these exciting times where discoveries of brown dwarfs are finally breaking 
through, model atmospheres are also rapidly becoming up to the task of inter- 
preting the observations and deriving new search strategies. Uniform grids of 
dwarf stars and brown dwarfs model atmospheres exist that extend from the tip 
to the toes of the main sequence - and beyond: 9000K to 900K, logg= 3.0-6.0, 
and [M/H]= 0.0 to —2.0 for the NextGen models. These large model grids al- 
lowed the construction of consistent interior and evolution models for VLMs that 
yield unprecedent agreement with globular cluster main sequences observed to 
0.1 Mq with HST. They led to the derivation of the important mass- luminosity 
relation for halo brown dwarfs and so to the realization that brown dwarfs cannot 
make up a significant fraction of the halo missing mass. 

The effective temperature scale of K to M type dwarfs with spectral types 
earlier than M6 is now unambiguously established, with only small uncertainties 
remaining from a possible incompleteness of existing TiO line lists. Grain forma- 
tion has been identified as an important process in M dwarfs and brown dwarfs 
atmospheres which could explain the long-standing difficulties of the models to 
reproduce the spectral distribution of dwarfs later than about M6. The results 
of the models indicate that it may not longer be assumed that the convection 
zone extends to the photosphere of late-type red dwarfs and brown dwarfs, and 
that their photospheric lithium abundance reflect their core temperature and 
mass. The basic assumption supporting the lithium test of substellarity is only 
valid for young, hot brown dwarfs such as those found in the Pleiades cluster. 
Fortunately, if the lithium test cannot identify transition objects and brown 
dwarfs of the field, the OS molecular opacity treatment and grain formation 
have introduce new gravity (hence age) effects in the NextGen models that were 
not seen in the previous Extended models and that will potentially allow to 
separate younger transitional objects from field stars as readily as from their 
location in color-color diagrams. For this the colors of late-type red dwarfs need 
to be known with good accuracy i.e. better than about 0.05 magnitude, which 
we find is not the case of many known late M dwarfs such as G1406, VBIO, and 
especially LHS2924. 

As cooler dwarfs are being discovered, spectral types are stretching far 
beyond the classical Morgan &; Keenan scheme. The lack of TiO bands in the 
optical, and the emergeance of CH4 opacities in the infrared in GD165B and 
G1229B call for an extension of the MK system beyond M9 to another spectral 
class (see D. Kirkpatrick, this volume). While the spectral class should only 
reflect the effective temperatures and not necessarily the mass of the objects, 
perhaps a suitable class for these objects would nevertheless be "T dwarfs" 
as in reminescence of J.C. Tarter who introduced the term "brown dwarf" now 
commonly accepted to designate substellar dwarfs, and Takashi Tsuji who led the 
field of late dwarfs atmospheres since the early 1960's, first introduced methane 
as a spectral indicator of substellarity, and who is retiring soon. Another spectral 
class, perhaps "P", will then be needed for dwarfs cooler then the condensation 
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point of water vapor including planets. In any case, studies of the optical spectra 
of Gliese 229B, GD165B, the DENIS and 2MASS objects and other late-type 
dwarfs will soon allow to determine the stellar surface coverage of dust clouds 
if such are present, and to verify if intrinsic spectral-type variability afflict cool 
dusty dwarfs. Models will be the subject of further investigations relative to 
grain formation and its effect on late-type dwarfs until they can reproduce the 
lower main sequence and lead the way into the regime of cool brown dwarfs. 
Finally, if brown dwarfs are not abundant in the halo, they certainly are in the 
galactic disk and their study remains one that shall flourish as the census of the 
solar neighborhood continues and the gap between planets and stars fills in. 
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